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Abstract 

The negative energy density of Casimir systems appears to violate general relativity energy 
conditions. However, one cannot test the averaged null energy condition (ANEC) using standard 
calculations for perfectly reflecting plates, because the null geodesic would have to pass through the 
plates, where the calculation breaks down. To avoid this problem, we compute the contribution to 
ANEC for a geodesic that passes through a hole in a single plate. We consider both Dirichlet and 
Neumann boundary conditions in two and three space dimensions. We use a Babinet's principle 
argument to reduce the problem to a complementary finite disk correction to the perfect mirror 
result, which we then compute using scattering theory in elliptical and spheroidal coordinates. In 
the Dirichlet case, we find that the positive correction due to the hole overwhelms the negative 
contribution of the infinite plate. In the Neumann case, where the infinite plate gives a positive 
contribution, the hole contribution is smaller in magnitude, so again ANEC is obeyed. These 
results can be extended to the case of two plates in the limits of large and small hole radii. This 
system thus provides another example of a situation where ANEC turns out to be obeyed when 
one might expect it to be violated. 
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I. INTRODUCTION 



The standard Casimir calculation of the energy density between a pair of parallel plates 
(see for example 1]) yields a negative energy density between the plates. While this result 
poses no problem in the calculation of the usual Casimir force, it presents a puzzle for 
general relativity. One can construct a spacetime with an arbitrary geometry Rxu simply by 
constructing the energy-momentum tensor to solve Einstein's equations 

T\u = -—^ ( Rxu — -^gxuR] ■ (1) 



8nG V 2^ 

The only way to prevent the appearance of exotic phenomena, such as closed timelike curves 
0, traversable wormholes or superluminal travel j^, is to place restrictions on the allowed 
energy-momentum tensors Tx,^. While these conditions are all obeyed in classical physics, 
the negative energy density of the quantum Casimir system violates most such conditions, 
including the weak energy condition (WEC) and the null energy condition (NEC), which 
require that T\yV'^V'^ > for timelike and null vectors respectively. A still weaker condition, 
which is still strong enough to rule out exotic phenomena (and to prove singularity theorems 
I^; IE 1^); IS that the null energy condition hold only when averaged over a complete 
geodesic (ANEC).^ Geodesies parallel to the plates obey NEC, so any candidate for ANEC 
violation would need to pass through the plates themselves. Therefore one cannot test ANEC 
using the standard Casimir calculation in which one imposes ideal boundary conditions, since 
this calculation is not valid within each plate. One approach to resolve this question is to 
model the plate using a domain wall background 0, [lO]; in this case the effect of the domain 
wall modifies the calculation significantly, so that ANEC is obeyed. In a number of other 
examples in which one might expect to find that ANEC is violated, explicit calculation shows 
that it is obeyed 11, 3|. Other calculations also show that energy condition violation is 
more difficult to achieve in realistic situations than idealized models would suggest 3, ^ ■ 
ANEC is also known to be obeyed by free scalar (15| and electromagnetic |[1Q] fields in flat 
spacetime. Other works have found restrictions on energy condition violation in flat space 



In this paper we consider an alternative modification to the Casimir problem that one 
might expect would allow the NEC violation between the plates to extend to ANEC viola- 
tion: we imagine a plate with a small hole, through which the geodesic can pass without 
encountering the material of the plate. Our primary calculation is for the case of a single 
Dirichlet plate, which also leads to a negative energy density for minimal coupling. We 
consider both two and three spatial dimensions. We also give extensions of this result to 
Neumann boundaries and to the case of two plates in extreme limits. 



II. NULL ENERGY CONDITION OUTSIDE A BOUNDARY 



Since our geodesic never passes through the material that actually imposes the boundary, 
we need only consider the quantum field in empty space. The effect of the boundary will 
be to modify the normal mode expansion for 0. We will then integrate V^VTxi,, where 
Txu is the stress-energy tensor, over the null geodesic perpendicular to the plate, passing 



For singularity theorems, the average is only over the future of the trapped surface. 
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through the center of the hole. The stress-energy tensor for a minimally-coupled scalar field 
is 

Ta. = 9a09,0 - ^vxu [d^<pdx(P] . (2) 

For a null vector, rjxuV'^V'^ = 0, so we have 

TxuV^V = . (3) 

For a static system, Tq^ = for i = 1,2, 3. If we further choose spatial coordinates in which 
Tij is diagonal and V = (1, v), then 

TxuV'V = 0' + 5^ ivA(l>f . (4) 

i 

Let the center of the hole lie at the origin and let the z axis be the direction perpendicular 
to the plate, along which the geodesic lies. For that path. 



III. BABINET'S PRINCIPLE 

It will greatly siraplify the scattering theory techniques we would like to use in our 
Casimir calculation 21, i^l to be able to consider a boundary condition in a local 
region. Therefore we apply a Babinet's principle argument to reexpress the result for a plate 
with a hole in terms of the results when the boundary condition is applied to an entire plate 
and when a complementary boundary condition is applied to a disk. The former is well- 
known, while the latter can be computed using scattering theory in elliptical or spheroidal 
coordinates. 

We start in empty space, and write the field there in terms of normal modes that are 
even or odd in the coordinate across the boundary. Next we consider a perfectly reflecting 
Dirichlet boundary with no holes. The free-space odd modes obey the boundary conditions, 
but the even modes do not. Instead, we have new even modes, which are just the odd mode 
on the right and minus the odd mode on the left, as shown in Fig. ^ If we let E denote a 
sum over the free-space even modes and O the same sum over the odd modes, then in free 
space we have E + 0, whereas with the barrier we have + 0. Therefore, the renormalized 
energy, the difference between the energy with the barrier and the energy in free space, is 
O — E. If we have Neumann conditions instead, the situation is precisely reversed and the 
energy is E — O. 

Now we consider the Dirichlet plate with holes of arbitrary shapes. Once again the odd 
modes are unaffected, and we have new even modes that vanish on the barrier but are 
continuous in the holes, as shown in Fig. |21 Since they are even, they satisfy Neumann 
conditions in the hole. Let us call the contribution of those modes A. The energy with the 
perforated barrier is thus A + O, so the renormalized energy is A — E. 

Finally, suppose that there are Neumann patches where the holes were. The even modes 
are unaffected, but there are new odd modes. In order to be odd and continuous they must 
satisfy Dirichlet conditions on the plane outside the patches. Thus, except for a change of 
sign on one side, these are the exact same modes of the previous paragraph, as shown in 
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FIG. 1: In free space (left) we have even and odd normal mode wave functions. In the presence 
of a Dirichlet plate (right), the even functions are replaced by the odd functions with a change of 
sign crossing the plate. 



FIG. 2: If there is a hole in the Dirichlet plate (left), the new even functions satisfy Neumann 
conditions in the hole and Dirichlet conditions elsewhere. The odd functions for a patch with the 
same shape as the hole (right) are the same except for a change of sign between sides. 

Fig. 121 Therefore the total energy is ^4 + and the renormalized energy is A — O. Thus we 
conclude that 

[Dirichlet plate w/hole] — [Complementary Neumann disk] = [Entire Dirichlet plane] (6) 
and similarly 

[Neumann plate w/hole] — [Complementary Dirichlet disk] = [Entire Neumann plane]. (7) 

IV. LINE SEGMENT IN TWO SPATIAL DIMENSIONS 

In this section we consider a scalar field in 2+1 dimensions with boundary conditions 
imposed on a line segment from x = —d to x = d. In circular coordinates, we can decompose 
a free, real, massless scalar field in modes as 




oo 
m=0 




k 
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where the prime on the summation sign indicates that for m = 0, there is no sin mode and 
instead of cosO = 1 we have 1/72. 

Now we go to elhptical coordinates. For notational consistency with the three-dimensional 
case, we consider the x-z plane. The foci will be located at x = id, and the the geodesic 
will run along the z axis. Elliptical coordinates /i, 9 are given by 

X = cosh /X cos 6* (9) 
z = (i sinh sin 6' (10) 

so 



r = = df-^^^^^^^±^ ^ ^e^ (11) 

as cxD. 

We define our Mathieu functions following the conventions of Abramowitz and Stegun 
(23! but extending their notation to be more similar to that of Bessel functions. The angular 
functions are cem{0,q) and sem{0,q). They satisfy 

?/" + (a-2gcos2^)?/ = (12) 

where q = {dk/2)'^, and are normalized so that 

/ deceUO,q)''= / deseUO,qf = n. (13) 
Jo Jo 

This normalization holds even for ceo, but there is no such function as scq. Thus the 
normalization is precisely the same as the circular functions used above, including the special 
case for m = 0. 

The radial functions of the first kind are Jem(/i, <?) and Jom(/i, g) and are precisely the 
functions called Mcm{fi,q) and Msm{fi,q) respectively in ji^]. They satisfy 

y" - (a - 2g cosh 2/i)y = (14) 

and go asymptotically to Jm{y/qe^) = Jm{kr). Note that the functions Jcm and Jom defined 



m 



2^ 1 have an additional factor of /2. Analogously, we denote the radial functions of 
the second kind as Yem{fJ',q) and Yom{fJ',q), which are McmifJ^yq) MsmifJ^yq) [2^ - 



The normalization of radial functions depends only on their asymptotics, so Jcm and 
Jom have the same normalization as Jm- Thus the field becomes 



m=0 

(15) 

where the prime on the summation sign indicates that the second term is included only for 
m > 0. 

Now we consider Neumann conditions along the line segment, which is /i = 0. The even 
wavefunctions obey the conditions already, because dJcm/dfi = at /i = 0, but the odd 
functions need to be modified. Instead of ^/'o^ = we need dipo^/dfj, = 0. For the free case 
we had 

Jorr^ifi, q) = l [Ho^^\fi, q) + Ho'g\^^, q)] (16) 
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where ifo*^^^ = Jom + iYom is the function called Msm in ji^] and Hom = Jom — iYom is 
Msm . Now we need 

i^oUf^, q) = \ [e'''Hog\f,, q) + Hog\^,, q)] (17) 



where 



iyo^^^'(0,g) 



and the derivative is with respect to /i. 

We can then compute the renormalized vacuum expectation value of the time-derivative 
term in Eq. (0), 

°^ poo 

(0^) = 7^y2 dkkuj {\ipOm{n,q)\'^ - Jom{n,qf) sem{n,qf (19) 

m=l 

and we can write 

|^o^(/i,g)P - Jo^(/i,g)^ = \ [{e''^' - l) //o«(^,g)2 + (e'^^^ - l) Hog\f^,qr] . (20) 

We want to extend the range of integration to include negative k. The term sem(/i, g)^ 
is unchanged by going to negative k, while Jo and Ho behave just like the corresponding 
Bessel functions. Thus the situation is exactly as in the circular case 21|: extending the 
range of integration exchanges the two terms in Eq. pO|l . so we can consider only the first 
term integrated over the entire real axis. When we close the contour at infinity, we get the 
contribution from the branch cut in c<j = v^- With k = in, the angular function becomes 
sCmiG, —q), and the radial functions have the same continuation as Bessel functions, 

Jom{fi,-q) = i"^IOm{fJ',q) (21) 
Hog\f^,-q) = -r(™+i)iro^(/x,g) (22) 

IT 



where lom and Kom are exactly as in f23j- Thus 



2i5^iiK) _ 1 _ / \m+l - -^Om(O^V^) /r,o\ 



where ip = ((iK/2)^ = —q. 

Putting it all together we have 



^ /»oo 

12^ = V / dK {e'^' - 1) K^Komif^, ^fseUO, -V) 



-2lir d^Pf^y^oUl^, ^fseUe, ■ (24) 



m=l 

On the axis, terms with m even vanish, so we have 
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where the prime on the summation sign indicates that wc sum over odd values of m. 

The other vacuum expectation value that we need is (((9^0)^). On the z axis, dzip is just 
the component of the gradient in the n direction, which differs from by the inverse of 
the metric coefficient 



cosh — cos 29 



h — dy cosh // — cos^ 6 — dy . (26) 

The calculation is otherwise similar. Instead of two powers of oo from time differentiation 
we just have the radial function differentiated with respect to /i, 

If instead we have Dirichlet conditions on the line segment, the odd functions will be un- 
modified, but for the even functions we need ipe^ = at = 0, so we have 

^eUl^, g) = ^ [e^^'i/e«(//, q) + Heg\ix, q)] (28) 



with 

Hei\Q,q) 



e- = -^^ (29) 



so 

Thus on the axis we have 

^^") --212* r d^P^.^'KeU^, -^f (31) 

^ 1^ Jo Ke^[0,ip) 

where the star on the summation sign indicates that we sum over even values of m. 



V. DISK IN THREE SPATIAL DIMENSIONS 

In this section we consider a scalar field with boundary conditions imposed on a disk of 
radius d in the x-y plane, centered at the origin. In spherical coordinates, we can decompose 
a free, real, massless scalar field in modes as 

OO I »oO I. 

0(r, ^, 0) = V V / dk -^ji{kr)Yimie, 0)e^'^* + c.c. (33) 
where ji is the spherical Bessel function. 
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Next we go to oblate spheroidal coordinates, given by 



X = d\/ {C,'^ + 1)(1 — 1]'^) cos (p = dcosh fi sin 9 cos (34) 
y = d\/ {^'^ + — 1]'^) sin (f) = d cosh ^ sin 9 sin (p (35) 
z = dr]^ = d sinh fi cos 9 (36) 

where (f) is the azimuthal angle, 1] = cos ^ is the coordinate akin to polar angle, and ^ = sinh n 
is the radial coordinate, with r = d^ for large ^. 

We define the prolate angular spheroidal function Sl^{c\ri) = Sn^^\c;ri) with c = kd, 
using the normalization of Meixner and Schafke ^J. Prolate spheroidal functions can be 
converted to the oblate ones appropriate to our situation hj k ^ ik and ^ —iC,. Thus 
our angular functions are S^{ic\ri), obeying the orthonormality relation 

' S-{^c■ r^)S^,{^c■ vfdv = iT^r^^.^nn' (37) 
I 2n + 1 [n — my. 

Using these functions we can define oblate spheroidal harmonics by analogy with spherical 
harmonics, 

y 47r (n + mj! 
obeying the analogous orthonormality relation 

' dr^ r d(Py::{lC- V, 0)*3^^'(^C; V, 0) = ^nn'Smm' ■ (39) 

-1 Jo 

We also define the radial spheroidal functions Rn^^\ic; —i^) and Rn^'^\ic; —i^), normal- 
ized by 

limi?™«(zc;-zO = Jn{c^) (40) 

5— >oo 

limR^^'\ic;-tO = Vnic^). (41) 

4— >oo 

The radial functions thus have the same normalization as spherical Bessel functions. Thus 
the field becomes 

0(e,^,0) = EE / dk^y!^{^c■,vA)K^'\^c■,-^Oe'-' + c.c. (42) 

If m + n is even, then S^{ic]r]) is an even function of r] and {d/d^)Rn^^\ic; —i^) = at 
^ = 0. Thus such wave functions will be continuous across the disk rj = 0. Similarly, if 
m + n is odd, then S^{ic\ rf) is an odd function of r] and R!^''^\ic\ 0) = 0, so the product is 
once again continuous. The i?^^-* functions do not have these boundary conditions, so they 
cannot be used in the vacuum wave functions. 

Now we consider Neumann conditions on the disk ^ = 0. If m + n is even, the functions 
obey the conditions already. Otherwise we need to combine i?^*^^'' and i^IT*^^'* to give the 
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desired condition. With i^JT^^^ = R^^^^ + iR^^'^^ and i?^^^^ = R^^^^ — iR^^^^ we can write 
the desired radial function 

-iO = \ [e'''^''^R^^'\ic; -i^ + R':^'\ic; -i^] (43) 

with the condition 

= (44) 



where the derivative is with respect to the second argument. 

The vacuum expectation value of the time derivative term, subtracting the free vacuum, 
then becomes 

n=0 m=—n ^ 

= rE E 7 dkk'u;\y:\tc;v,<P)\'[{e'^'(-^ -1) Rf'\tc;-tO' 

n=0 m=-n 

+ (e-'^^(^^) - 1) (ic; -iO'] (45) 



where the prime on the summation sign means that only odd values of m + n are included. 

We want to extend the range of integration to include negative k, which changes the sign 
of c. We can implement this change by changing the sign of ^ in Rn''^^ and Rn^'^\ since 
these real functions are not affected by complex conjugation. The functions i?™''^'' and i^JT*'^'* 
have opposite parity under this transformation, so i?™''^'' and i?™*'"''' change places. Thus 
including negative k in the first term gives the second term, and we have 



^ oo n 

('^') = rE E7 dkkMy':{ic-,vA)?{e'''^-^ -i)R':^'\ic-,-ie ■ m 

™_ „ J — oo 



n=0 m=—n 



If we take k and thus c in the upper half plane we will get spheroidal functions whose 
parameter goes to negative real infinity. Therefore we can close the contour at infinity and 
obtain an integral along the branch cut on the imaginary axis associated with the square 
root in a;. 



- ^E E 7 dKK^y-{r,v,<P)ne'''^^^ - I) R-^'\r,-^0' 

nrt A — "0 



n=0 m=—n 



^ oo n „oo n^l 



n=0 m=—n 

where 'j = ic = ikd = —nd. On the axis, we have 

?o(i) 



('^') = -^E' r^^^' S(3)!!^'!!! iyn(7;i,0)pi?fH7;-^O^ (48) 

where we have specialized to m = because the contributions from nonzero m vanish on 
the axis, leaving only a sum over odd values of n. Similarly, on the axis we have 

ii^^^r) = -^E' r^--|M!|^l>^^(7;i,0) (49) 
TTfl Jo R^: "{r, 0) 
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where the primes on the radial functions indicate derivatives with respect to the second 
argument, and the metric coefficient 



5x 



d 



^2 + ^2 



(50) 



becomes equal to d on the axis. 

For Dirichlet conditions on the disk, we have the analogous results 




(51) 



(52) 



where we have again specialized to the axis so that only m = contributes, the derivatives 
of the radial functions are again with respect to the second argument, and the star on the 
summation sign indicates that now we sum over even values of n. 

VI. NUMERICAL CALCULATION 

For a null geodesic V^'^ perpendicular to the plate and passing through the center of the 
hole, the contribution to ANEC is given by Eq. (0). We can compute the results for the 
complementary disk using the formulae derived in the previous sections. We then add the 
complete plate results. 



with Dirichlet conditions giving the upper sign and Neumann the lower. 

In two dimensions, we compute the Mathieu functions using the package of Alhargan 
(iol I27I I , with some minor modifications: we use 80-bit double precision throughout the cal- 
culation to accommodate the extreme dynamical range needed for the wide range of Mathieu 
function parameters we use, and we have adapted the code to use our set of normalization 
conventions. These C++ routines are then imported into Mathematica, where we can use 
efficient routines for numerical sums and integrals. 

In three dimensions, we use the Mathematica spheroidal harmonic package of Falloon ji^ . 
We have updated it to fix incompatibilities with the latest version of Mathematica and to 
avoid memory leaks. We have also made a number of efficiency optimizations appropriate to 
the unusual demands we make on the code (for example, we changed the caching structure 
so that it is appropriate to the way we call the functions, with the same arguments but 
different parameters rather than the other way around; we also wrote specific code for the 
modified radial function of the third kind to avoid cancellations of exponentially growing 
quantities). 

Figure El shows the contributions to NEC for Dirichlet plates with holes of unit radius in 
two and three spatial dimensions, as functions of distance along the axis. Using Babinet's 




in three dimensions 



in two dimensions, and 



(53) 
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FIG. 3: Contributions to NEC in two dimensions (left) and three dimensions (right) for a Dirichlet 
plate with a hole of unit radius, as functions of distance along the axis passing through the center 
of the hole. Extrapolation is used for points at a distance less than 0.15 in the left panel and 0.25 
in the right panel. The dotted lines show the perfect mirror result. 



principle, we have computed the sum of contributions from a Neumann disk and an infinite 
Dirichlet mirror. At small distances, the contributions from both the finite disk and the 
infinite mirror diverge like 1/z""^^, where z is the distance from the origin and n is the 
spatial dimension. The true result, however, does not diverge (the origin is just a point in 
empty space) and by symmetry must have zero slope at the origin. This cancellation provides 
a highly nontrivial check on our calculation. Going all the way to the origin would require 
infinite precision; in two dimensions we stop at a distance O.lSd, while in three dimensions 
we stop at distance 0.25(i. At these values, our curves already show this cancellation clearly. 
We also extrapolate our result (without putting in any restrictions on the extrapolation at 
the origin) and find that it goes smoothly to a finite value with zero slope at 2; = 0. A less 
stringent check is that our calculation approaches the perfect mirror at large distances; in 
our approach this result simply tells us that the finite disk contribution is going to zero fast 
enough. 



VII. DISCUSSION 

The results in Fig. |21 are quite striking. In both cases, far from the origin we see the 
negative contribution to ANEC that we would expect from the standard calculation. Near 
the plate, however, the hole leads to a large positive contribution. Integrating the results 
shown in Fig. El gives a total contribution (including both sides of the plate) of 1.63 x 
10~^/(i^ in three dimensions and 4.53 x 10~^/(i^ in two dimensions. The positive contribution 
overwhelms the negative contribution, so that ANEC is obeyed. 

We can look at these results from a different point of view by considering conformal 
coupling. Since the NEC contribution in this case differs from that in minimal coupling by 
a total derivative, it leads to the same results for ANEC. From this point of view, one might 
not expect any ANEC violation in the case at hand, because the quantum contributions to 
the perfect mirror vanish. (In the perfect mirror case, changing from minimal to conformal 
coupling effectively moves the negative contributions from the region outside the boundary 
onto the boundary itself jl^.) However, from this point of view one might just as well expect 
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FIG. 4: Contributions to NEC in two dimensions (left) and three dimensions (right) for a Neumann 
plate with a hole of unit radius, as functions of distance along the axis passing through the center 
of the hole. Extrapolation is used for points at a distance less than 0.11 in the left panel and 0.25 
in the right panel. The dotted lines show the perfect mirror results. 

Neumann conditions to violate ANEC; while the perfect Neumann mirror result is positive 
in minimal coupling, it also vanishes for conformal coupling. The results for Neumann 
conditions are shown in Fig. EJ Once again integration gives positive results, 2.10 x 10~^/d^ 
in three dimensions and 1.77 x lO^^/rf^ in two dimensions, so ANEC is obeyed. Thus again 
from this point of view, one finds ANEC obeyed more often than would be naively expected. 

A conformal field between two plates would have a constant negative energy density. We 
can estimate how our results would extend to the case of parallel plates with holes, in two 
limits. If the separation between the plates i is much smaller than the radius of the hole d, 
the two plates are equivalent to a single plate and ANEC continues to hold. In the other 
extreme, if the separation between the plates is large compared to the radius of the hole, 
then we can assume that the change A in ANEC induced by adding a hole in one plate is 
unaffected by the other plate. In three dimensions, we obtain a contribution to ANEC of 



where the first term is the effect of the each of the two holes individually and the second 
term is the standard contribution from the two plates. For Dirichlet conditions, A = 1.63 x 
10~^/d^. Thus Eq. (jSD) gives a positive result as long as £ > 1.6d, which surely includes its 
entire range of applicability. Similar results hold for the other cases we have considered. 
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